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CD Abstract. The radio millisecond pulsar PSR B1620— 26 is part of an 

extraordinary triple star system in the globular cluster M4. The inner 
companion to the neutron star is thought to be a white dwarf of mass 



o 

(N 



qj | mi — 0.3 Mq in an orbit of period ~ 0.5 yr. The nature and orbital 

/~\ . characteristics of the second, more distant companion, have remained a 

mystery for many years. A theoretical analysis of the latest available 
radio pulsar timing data is presented here, allowing us to determine ap- 
proximately the mass and orbital parameters of the second companion. 
j> ■ Remarkably, the current best-fit parameters correspond to a second com- 

panion of planetary mass, with m2sini2 — 7 x 1O _3 M0, in an orbit of 
■ eccentricity e2 — 0.45 and with a large semimajor axis 02 — 60 AU. The 

short dynamical lifetime of this very wide triple in M4 suggests that large 
numbers of such planets must be present in globular clusters. We also 
£^ I address the question of the anomalously high eccentricity of the inner 

binary pulsar. While this eccentricity could have been induced during 
the same dynamical interaction that created the triple, we find that it 
Qh' could also naturally arise from long-term secular perturbation effects in 

q the triple, combining the general relativistic precession of the inner orbit 

$_i ■ with the Newtonian gravitational perturbation by the outer planet. 
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Introduction 



PSR B1620— 26 is a unique millisecond radio pulsar. The pulsar is a member 
of a hierarchical triple system located in or near the core of the globular cluster 
M4. It is the only radio pulsar known in a triple system, and the only triple 
system known in any globular cluster. The inner binary of the triple contains 
the ~ 1.4 M neutron star with a ~ 0.3 Mq white-dwarf companion in a 191- 
day orbit (Lyne et al. 1988; McKenna & Lyne 1988). The triple nature of the 
system was first proposed by Backer (1993) in order to explain the unusually high 
residual second and third pulse frequency derivatives left over after subtracting 
a standard Keplerian model for the pulsar binary. 

The pulsar has now been timed for 12 years since its discovery (Thorsett, Ar- 
zoumanian, & Taylor 1993; Backer, Foster, & Sallmen 1993; Backer & Thorsett 
1995; Arzoumanian et al. 1996; Thorsett et al. 1999). These observations have 
not only confirmed the triple nature of the system, but they have also provided 
constraints on the mass and orbital parameters of the second companion. Earlier 
calculations using three pulse frequency derivatives suggested that the mass of 
the second companion could be anywhere between ~ 1CP 3 — 1-M©, with corre- 
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sponding orbital periods in the range ~ 10 2 - 10 3 yr (Michel 1994; Rasio 1994; 
Sigurdsson 1995). More recent calculations using four frequency derivatives and 
preliminary measurements of the orbital perturbations of the inner binary fur- 
ther constrained the mass of the second companion, and suggested that it is most 
likely a giant planet or a brown dwarf of mass ~ 0.01 M at a distance ~ 50 AU 
from the pulsar binary (Arzoumanian et al. 1996; Joshi & Rasio 1997, here- 
after JR97). In §2 below we summarize our results from the latest theoretical 
analysis of the pulsar timing data (Ford et al. 2000; hereafter FJRZ), including 
the most recent observations of Thorsett et al. (1999; hereafter TACL). The 
data now include measurements of five pulse frequency derivatives, as well as 
improved measurements and constraints on various orbital perturbation effects 
in the triple. 

Previous optical observations by Bailyn et al. (1994) and Shearer et al. 
(1996) using ground-based images of M4 had identified a possible optical coun- 
terpart for the pulsar, consistent with a ~ 0.5 M Q main-sequence star, thus 
contradicting the theoretical results, which suggest a much lower-mass com- 
panion. However, it also seemed possible that the object could be a blend of 
unresolved fainter stars, if not a chance superpositon. Later HST WFPC2 obser- 
vations of the same region by Bailyn (private communication) have resolved the 
uncertainty. The much higher resolution (~ 0.1") HST image shows no optical 
counterpart at the pulsar position, down to a magnitude of V ~ 23, therefore 
eliminating the presence of any main-sequence star in the system. 

PSR B1620— 26 is not the first millisecond pulsar system in which a planet 
(or brown dwarf) has been detected. The first one, PSR B1257+12, is a single 
millisecond pulsar with three clearly detected inner planets (all within 1 AU) of 
terrestrial masses in circular orbits around the neutron star (Wolszczan 1994). 
Preliminary evidence for at least one giant planet orbiting at a much larger dis- 
tance from the neutron star has also been reported (Wolszczan 1996; JR97). 
In PSR B1257+12 it is likely that the planets were formed in orbit around the 
neutron star, perhaps out of a disk of debris left behind following the complete 
evaporation of a previous binary companion (see, e.g., Podsiadlowski 1995). 
Such an evaporation process has been observed in several eclipsing binary mil- 
lisecond pulsars (e.g., Nice et al. 2000), where the companion masses have been 
reduced to ~ 0.01 M & by ablation. These companions used to be ordinary white 
dwarfs and, although their masses are now quite low, they cannot be properly 
called either planets or brown dwarfs (Rasio, Pfahl, & Rappaport 2000). 

In PSR B1620— 26, the hierarchical triple configuration of the system and 
its location near the core of a dense globular cluster suggest that the second 
companion was acquired by the pulsar following a dynamical interaction with 
another object (Rasio, McMillan & Hut 1995; Sigurdsson 1993, 1995; FJRZ). 
This object could have been a primordial binary with a low-mass brown-dwarf 
component, or a main-sequence star with a planetary system containing at least 
one massive giant planet. Indeed the possibility of detecting "scavenged" plan- 
ets around millisecond pulsars in globular clusters was discussed by Sigurdsson 
(1992) even before the triple nature of PSR B1620— 26 was discovered. Several 
versions of such a dynamical formation scenario are possible for the triple sys- 
tem, all involving dynamical exchange interactions between binaries in the core 
of M4. In the most likely scenario, studied in detail by FJRZ, a pre-existing 
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binary millisecond pulsar has a dynamical interaction with a wide star-planet 
system, which leaves the planet bound to the binary pulsar while the star is 
ejected. From numerical scattering experiments we found that the probability 
of retaining the planet, although smaller than the probability of retaining the 
star, is always significant, with a branching ratio ~ 10% — 30% for encounters 
with pericenter distances r p /a,i ~ 0.2 — 1, where ctj ~ 50 AU is the typical ini- 
tial star-planet separation. All the observed parameters of the triple system 
are consistent with such a formation scenario, which also allows the age of the 
millisecond pulsar (most likely ^ 10 9 yr) to be much larger than the lifetime of 
the triple (as short as ~ 10 7 yr if it resides in the core of the cluster). 

Objects with masses ~ 0.001 — 0.01 M & have recently been detected around 
many nearby solar-like stars in Doppler searches for extrasolar planets (see 
Marcy & Butler 1998 for a review). In several cases (e.g., v And), more than 
one object have been detected in the same system, clearly establishing that they 
are members of a planetary system rather than a very low-mass stellar (brown 
dwarf) binary companion. For the second companion of PSR B1620— 26, of mass 
~ 0.01 Mq, current observations and theoretical modeling do not make it possi- 
ble to determine whether the object was originally formed as part of a planetary 
system, or as a brown dwarf. In this paper, we will simply follow our prejudice, 
and henceforth will refer to the object as "the planet." 

One aspect of the system that remains unexplained, and can perhaps pro- 
vide constraints on its formation and dynamical evolution, is the unusually high 
eccentricity e\ = 0.025 of the inner binary. This is much larger than one would 
expect for a binary millisecond pulsar formed through the standard process of 
pulsar recycling by accretion from a red-giant companion. During the mass 
accretion phase, tidal circularization of the orbit through turbulent viscous dis- 
sipation in the red-giant envelope should have brought the eccentricity down to 
£ 10" 4 (Phinney 1992). At the same time, however, the measured value may 
appear too small for a dynamically induced eccentricity. Indeed, for an initially 
circular binary, the eccentricity induced by a dynamical interaction with another 
star is an extremely steep function of the distance of closest approach (Rasio 
& Heggie 1995). Therefore a "typical" interaction would be expected either 
to leave the eccentricity unchanged, or to increase it to a value of order unity 
(including the possibility of disrupting the binary) . 

Secular perturbations in the triple system can also lead to an increase in the 
eccentricity of the inner binary. A previous analysis assuming nearly coplanar 
orbits suggested that, starting from a circular inner orbit, an eccentricity as large 
as 0.025 could only be induced by the perturbation from a stellar-mass second 
companion (Rasio 1994), which is now ruled out. For large relative inclinations, 
however, it is known that the eccentricity perturbations can in principle be 
considerably larger (Kozai 1962; see Ford, Kozinsky & Rasio 2000, hereafter 
FKR, for a recent treatment). In addition, the Newtonian secular perturbations 
due to the tidal field of the second companion can combine nonlinearly with 
other perturbation effects, such as the general relativistic precession of the inner 
orbit, to produce enhanced eccentricity perturbations (see FKR and §3 below). 
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2.1. Pulse Frequency Derivatives 

The standard method of fitting a Keplerian orbit to timing residuals cannot 
be used when the pulsar timing data cover only a small fraction of the orbital 
period (but see TACL for an attempt at fitting two Keplerian orbits to the 
PSR B1620-26 timing data). For PSR B1620-26, the duration of the obser- 
vations (~ 10 yr) is short compared to the likely orbital period of the second 
companion, which is ^ 100 yr. In this case, pulse frequency derivatives (coeffi- 
cients in a Taylor expansion of the pulse frequency around a reference epoch) 
can be derived to characterize the shape of the timing residuals (after subtrac- 
tion of a Keplerian model for the inner binary) . It is easy to show that from five 
well-measured and dynamically-induced frequency derivatives one can obtain a 
complete solution for the orbital parameters and mass of the companion, up to 
the usual inclination factor (JR97). 

JR97 used the first four time derivatives of the pulse frequency to solve 
for a one-parameter family of solutions for the orbital parameters and mass 
of the second companion. The detection of the fourth derivative, which was 
marginal at the time, has now been confirmed (TACL). In addition, we now 
also have a preliminary measurement of the fifth derivative. This allows us in 
principle to obtain a unique solution, but the measurement uncertainty on the 
fifth derivative is very large, giving us correspondingly large uncertainties on 
the theoretically derived parameters of the system. Equations and details on 
the method of solution are presented in JR97, and will not be repeated here. 

Our new solution is based on the latest available values of the pulse fre- 
quency derivatives, obtained by TACL for the epoch MJD 2448725.5: 

Spin Period P = 11.0757509142025(18) ms 

Spin frequency/ = 90.287332005426(14) s" 1 

/ = -5.4693(3) x 10~ 15 s- 2 

/ = 1.9283(14) x 10~ 23 s~ 3 

/ = 6.39(25) x 10" 33 s" 4 

= -2.1(2) x HT 4( V 5 

/( 5 ) = 3(3) x KT 49 s- 6 

Here the number in parenthesis is a conservative estimate of the formal la error 
on the measured best-fit value, taking into account the correlations between 
parameters (see TACL for details). It should be noted that the best-fit value 
for the fourth derivative quoted earlier by Arzoumanian et al. (1996) and used 
in JR97, /( 4 ) = —2.1(6) x 10 -40 s -5 , has not changed, while the estimated la 
error has decreased by a factor of three. This gives us confidence that the new 
measurement of f^ 5 \ although preliminary, will not change significantly over the 
next few years as more timing data become available. 

Since the orbital period of the second companion is much longer than that 
of the inner binary, we treat the inner binary as a single object. Keeping the 
same notation as in JR97, we let m\ = ttins + m c be the mass of the inner 
binary pulsar, with ttt-ns the mass of the neutron star and m c the mass of the 
(inner) companion, and we denote by m,2 the mass of the second companion. The 
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orbital parameters are the longitude A2 at epoch (measured from pericenter), the 
longitude of pericenter u>2 (measured from the ascending node), the eccentricity 
e2, semimajor axis 02, and inclination 12 (such that sin 12 = 1 for an orbit seen 
edge-on). They all refer to the orbit of the second companion with respect 
to the center of mass of the system (the entire triple). A subscript 1 for the 
orbital elements refers to the orbit of the inner binary. We assume that ttt-ns = 
1.35 Mq, giving m c sin i\ = 0.3 Mq, whereii is the inclination of the inner binary 
(Thorsett et al. 1993), and we take sinii = 1 for the analysis presented in this 
section since our results depend only very weakly on the inner companion mass. 

The observed value of / is in general determined by a combination of the in- 
trinsic spin-down of the pulsar and the acceleration due to the second companion. 
However, in this case, the observed value of / has changed from —8.1 x 10~ 15 s~ 2 
to —5.4 x 10~ 15 s~ 2 over 11 years (TACL). Since the intrinsic spin-down rate is 
essentially constant, this large observed rate of change indicates that the ob- 
served / is almost entirely acceleration-induced. Similarly, the observed value 
of / is at least an order of magnitude larger than the estimate of / from in- 
trinsic timing noise, which is usually not measurable for old millisecond pulsars 
(see Arzoumanian et al. 1994, TACL, and JR97). Intrinsic contributions to the 
higher derivatives should also be completely negligible for millisecond pulsars. 
Hence, in our analysis, we assume that all observed frequency derivatives are 
dynamically induced, reflecting the presence of the second companion. 

Figure 1 illustrates our latest one-parameter family of solutions, obtained 
using the updated values of the first four pulse frequency derivatives. There 
are no significant differences compared to the solution obtained previously in 
JR97. The vertical solid line indicates the unique solution obtained by including 
the fifth derivative. It corresponds to a second companion mass mi sin 12 = 
7.0 x 1O~ 3 M0, eccentricity e2 = 0.45, and semimajor axis 02 = 57 AU. For 
a total system mass m\ + 7772 = 1.65 M this gives an outer orbital period 
P 2 = 308yr. 

It is extremely reassuring to see that the new measurement of is con- 
sistent with the family of solutions obtained previously on the basis of the first 
four derivatives. The implication is that the signs and magnitudes of these five 
independently measured quantities are all consistent with the basic interpreta- 
tion of the data in terms of a second companion orbiting the inner binary pulsar 
in a Keplerian orbit. For comparison, the two vertical dashed lines in Figure 1 
indicate the change in the solution obtained by decreasing the value of by 
a factor of 1.5 (right), or increasing it by a factor of 1.5 (left). Note that lower 
values of give higher values for mi- If we vary the value of within its 
entire la error bar, all solutions are allowed, except for the extremely low-mass 
solutions with 7772 sin 12 <J 0.002 Mq. In particular, the present la error on 
does not strictly rule out a hyperbolic orbit (e.2 > 1) for 1712- However, it is still 
possible to derive a strict upper limit on 7712 by considering hyperbolic solutions 
and requiring that the relative velocity at infinity of the perturber be less than 
the escape speed from the cluster core. Indeed, for 777-2 sin 12 > 0.055 Mq, FJRZ 
find that the relative velocity at infinity would exceed the escape speed from the 
cluster core (~ 12kms _1 ; see Peterson et al. 1995). A strict lower limit on the 
mass, 7772 ^ 2 x 10~ 4 M©, can also be set by requiring that the orbital period 
of the second companion be longer than the duration of the timing observations 
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Figure 1. Allowed values of the semimajor axis <72, mass 777.2, eccentric- 
ity e2, longitude at epoch A2, and longitude of pericenter u>2 for the second 
companion of PSR B1620— 26, using the latest available values for four pulse 
frequency derivatives (left). On the right, we show the number of accepted 
realizations (N) of the triple for different values of 777,2 and the corresponding 
distance r\i of the second companion from the inner binary in our Monte Carlo 
simulations. Accepted realizations are those leading to short-term orbital per- 
turbation effects consistent with the current observations. The vertical solid 
lines indicates the complete solution obtained by including the preliminary 
value of the fifth derivative. The two dashed lines indicate the solutions ob- 
tained by decreasing or increasing the value of /( 5 ) by a factor of 1.5. 



(about 10 yr). Note that all solutions then give dynamically stable triples, even 
at the low-mass, short-period limit (Eggleton & Kiseleva 1995; JR97). 

2.2. Orbital Perturbations 

Additional constraints and consistency checks on the model can be obtained by 
considering the perturbations of the orbital elements of the inner binary caused 
by the presence of the second companion. These include a precession of the peri- 
center, as well as short-term linear drifts in the inclination and eccentricity. The 
drift in inclination can be detected through a change in the projected semimajor 
axis of the pulsar. The semimajor axis itself is not expected to be perturbed 
significantly by a low- mass second companion (Sigurdsson 1995). 

The latest measurements, obtained by adding a linear drift term to each 
orbital element in the Keplerian fit for the inner binary (TACL) give: 

Cj x = -5(8) x 10- 5 degyr -1 , 
ei = 0.2(1.1) x 10~ 15 s _1 , 
x p = -6.7(3) x IO" 13 , 
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where x p = a p smii is the projected semimajor axis of the pulsar. Note that 
only x p is clearly detected, while the other two measurements only provide upper 
limits. 

We have used these measurements in FJRZ to constrain the system by 
requiring that all our solutions be consistent with these secular perturbations. 
Specifically, we perform Monte Carlo simulations, constructing a large number of 
random realizations of the triple system in 3D, and accepting or rejecting them 
on the basis of compatibility with the measured orbital perturbations (see JR97 
for details) . The eccentricity of the outer orbit is selected randomly assuming a 
thermal distribution, and the other orbital parameters are then calculated from 
the standard solutions obtained in §2.1. The unknown inclination angles i\ and 
i2 are generated assuming random orientations of the orbital planes, and the two 
position angles of the second companion are determined using ii, 12, 0J2, A2, and 
an additional undetermined angle a, which (along with i\ and 12) describes the 
relative orientation of the two orbital planes. The perturbations are calculated 
analytically for each realization of the system, assuming a fixed position of the 
second companion. The perturbation equations are given in Rasio (1994) and 
JR97. 

Figure 1 (right panel) shows the resulting probability distributions for the 
mass m,2 and the current separation (at epoch) r\2 of the second companion. The 
Monte Carlo trials were performed using only our standard solution based on 
four pulse frequency derivatives, since the fifth derivative is still only marginally 
detected. The solid lines indicates the values given by the preliminary measure- 
ment of and assuming sin 22 = 1- The most probable value of m.2 — 0.01 M & 
is consistent with the range of values obtained from the complete solution us- 
ing the fifth derivative. We have also calculated the probability distribution 
of /( 5 ) predicted by our Monte Carlo simulations, and found that it is con- 
sistent with the preliminary measurement, providing yet another independent 
self-consistency check on our model. 

3. Secular Eccentricity Perturbations 

We now examine the possibility that the inner binary eccentricity was induced 
by the secular gravitational perturbation of the second companion. In the hi- 
erarchical three-body problem, analytic expressions for the maximum induced 
eccentricity and the period of long-term eccentricity oscillations are available in 
certain regimes, depending on the eccentricities and relative inclination of the 
orbits. 

3.1. Planetary Regime 

For orbits with small eccentricities and a small relative inclination (i ^ 40°), the 
classical solution for the long-term secular evolution of eccentricities and longi- 
tudes of pericenters can be written in terms of an eigenvalue and eigenvector 
formulation (e.g., Murray & Dermott 2000; see Rasio 1995 for simplified expres- 
sions in various limits). This classical solution is valid to all orders in the ratio of 
semimajor axes. The eccentricities oscillate as angular momentum is transferred 
between the two orbits. The precession of the orbits (libration or circulation) is 
coupled to the eccentricity oscillations, but the relative inclination remains ap- 
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proximately constant. In this regime, it can be shown that a stellar-mass second 
companion would be necessary to induce the observed eccentricity in the inner 
binary (Rasio 1994; 1995). Such a large mass for the third body has been ruled 
out by recent timing data (see §2), implying that secular perturbations from a 
third body in a nearly coplanar and circular orbit does not explain the observed 
inner eccentricity. 

Sigurdsson (1995) has suggested that it may be possible for the secular 
perturbations to grow further because of random distant interactions of the 
triple with other cluster stars at distances ~ 100 AU, which would perturb the 
long-term phase relation between the inner and outer orbits, allowing the inner 
eccentricity to "random walk" up to a much larger value. However, the current 
most probable solution from the timing data requires the second companion to 
be in a very wide orbit (separation ^ 40 AU), giving it a very short lifetime in 
the cluster, and leaving it extremely vulnerable to disruption by such repeated 
weak encounters (see JR97 for a more detailed discussion). 

3.2. High-Inclination Regime 

For a triple system formed through a dynamical interaction, there is no reason 
to assume that the relative inclination of the two orbits should be small. When 
the relative inclination of the two orbital planes is ^ 40°, a different regime 
of secular perturbations is encountered. This regime has been studied in the 
past using the quadrupole approximation (Kozai 1962; see Holman, Touma, & 
Tremaine 1997 for a recent discussion). Here the relative inclination i of the 
two orbits and the inner eccentricity e\ are coupled by the integral of motion 
6 = (1 — ef)cos 2 z (Kozai's integral). Thus, the amplitude of the eccentricity 
oscillations is determined by the relative inclination. It can be shown that large- 
amplitude eccentricity oscillations are possible only when B < 3/5 (Holman 
et al. 1997). For an initial eccentricity e\ ~ and initial inclination iq this 
implies i$ > cos -1 \/3/5 ~ 40° and the maximum eccentricity is then given 
by ei ma x = [1 — (5/3) cos 2 io} 1 ^ 2 , which approaches unity for relative inclinations 
approaching 90°. For a sufficiently large relative inclination, this suggests that it 
should always be possible to induce an arbitrarily large eccentricity in the inner 
binary, and that this could provide an explanation for the anomalously high 
eccentricity of the binary pulsar in the PSR B1620— 26 system (Rasio, Ford, 
Sz Kozinsky 1997). However, there are two additional conditions that must be 
satisfied for this explanation to hold. 

First, the timescale for reaching a high eccentricity must be shorter than 
the lifetime of the system. Although the masses, initial eccentricities, and ratio 
of semimajor axes do not affect the maximum inner eccentricity, they do affect 
the period of the eccentricity oscillations (see Holman et al. 1997; Mazeh et al. 
1996). The inner longitude of periastron precesses with this period, which can 
be quite long, sometimes exceeding the lifetime of the system in the cluster. The 
masses also affect the period, but they decrease the amplitude of the eccentric- 
ity oscillations only when the mass ratio of the inner binary approaches unity 
(FKR) . In Figure 2 we compare the period of the eccentricity oscillations to the 
lifetime of the triple in M4. It is clear that for most solutions the timescale to 
reach a large eccentricity exceeds the lifetime of the triple. The only possible 
exceptions are for very low-mass planets (rr»2 <^ 0.002 M Q ) and with the triple 
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residing far outside the cluster core. These cannot be excluded, but are certainly 
not favored by the observations (see §2). 

The second problem is that other sources of perturbations may become sig- 
nificant over these long timescales. In particular, for an inner binary containing 
compact objects, general relativistic effects can become important. This turns 
out to play a crucial role for the PSR B1620— 26 system, and we address this 
question in detail in the next section. 

3.3. General Relativistic Effects 

Additional perturbations that affect the longitude of periastron can indirectly 
change the evolution of the eccentricity of the inner binary in a hierarchical 
triple. For example, tidally- or rotationally-induced quadrupolar distortions, as 
well as general relativity, can cause a significant precession of the inner orbit 
for a sufficiently compact binary. If this additional precession is much slower 
than the precession due to the secular perturbations, then the eccentricity os- 
cillations are not significantly affected. However, if the additional precession is 
faster than the secular perturbations, then eccentricity oscillations are severely 
damped (see, e.g., Holman et al. 1997). In addition, if the two precession periods 
are comparable, then a type of resonance can occur that leads to a significant 
increase in the eccentricity perturbation (FKR). 

Figure 2 compares the various precession periods for PSR B1620-26 as a 
function of the second companion mass ni2 for the entire one-parameter family 
of standard solutions constructed in §2. The precession period Pmgh i for high- 
inclinations was calculated using the approximate analytic expression given by 
Holman et al. (1997, eq. 3), while -Plow i is from the classical solution (§3.1). The 
general relativistic precession period for the inner binary, Pgr, is also shown. 
Note that, although we have labeled the plot assuming sin 22 = 1, only .Plow i 
has an explicit dependence on mi (and it is calculated here for sini2 = 1). 

For nearly all solutions we find that the general relativistic precession is 
faster than the precession due to the Newtonian secular perturbations. The 
only exceptions are for low-mass second companions [ra^ <J 0.005 M©) in low- 
inclination orbits. However, we have already mentioned above (§3.1) that in this 
case the maximum induced eccentricity could not reach the present observed 
value. Most remarkably, however, we also see from Figure 2 that, for the most 
probable solution (based on the current measured value of and indicated by 
the vertical solid line in the figure), the two precession periods are nearly equal 
for a low-inclination system. This suggests that resonant effects may play an 
important role in this system, a possibility that we have explored in detail using 
numerical integrations. We describe the numerical results in the next section. 

3.4. Numerical Integrations of the Secular Evolution Equations 

If the quadrupole approximation used for the high-inclination regime is ex- 
tended to octupole order, the resulting secular perturbation equations approx- 
imate very well the long-term dynamical evolution of hierarchical triple sys- 
tems for a wide range of masses, eccentricities, and inclinations (FKR; Kry- 
molowski &: Mazeh 1999). We have used the octupole-level secular perturbation 
equations derived by FKR to study the long-term eccentricity evolution of the 
PSR B1620— 26 triple. We integrate the equations using the variables e\ sinwi, 
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Figure 2. Comparison of the various secular precession timescales in the 
PSR B1620— 26 triple, as a function of the mass of the second companion in 
the standard solution of §2 (left). Also shown for comparison is the lifetime 
of the triple, both in the cluster core (r c ) and at the half-mass radius (Th m )- 
On the right, we show the maximum eccentricity of the binary pulsar induced 
by secular perturbations in the triple (including the Newtonian perturbation 
by the second companion, and the general relativistic precession of the inner 
orbit), as a function of the mass of the second companion in the one-parameter 
family of solutions from §2. Note the clear resonant interaction between the 
Newtonian and relativistic perturbations. 

ei costal, e2 sina^, and &2 cos u>2, where u\ and u>2 are the longitudes of pericen- 
ters. This avoids numerical problems for nearly circular orbits, and also allows 
us to incorporate easily the first-order post-Newtonian correction into our in- 
tegrator, which is based on the Burlish-Stoer method. We assume that the 
present inner eccentricity is due entirely to the secular perturbations and that 
the initial eccentricity of the binary pulsar was much smaller than its present 
value. In addition, we restrict our attention to the standard one-parameter fam- 
ily of solutions constructed in §2.1. From the numerical integrations we can then 
determine the maximum induced eccentricity. 

In Figure 2 (right panel) we show this maximum induced eccentricity in the 
inner orbit as a function of the second companion mass for several inclinations. 
For most inclinations and masses, we see that the maximum induced eccentricity 
is significantly smaller than the observed value, as expected from the discussion 
of §3.3. However, for a small but significant range of masses near the most 
probable value (approximately 0.0055 M Q < rri2 < 0.0065 Mq), the induced 
eccentricity for low-inclination systems can reach values ^ 0.02. 

As already pointed out in §3.2, the maximum induced eccentricity may also 
be limited by the lifetime of the triple system. From Figure 2 we see that, near 
resonance, we expect Plow i — -Pgr — 10 7 yr, which is comparable to the lifetime 
of the triple in the cluster core, and much shorter than the lifetime outside 
of the core. For solutions near a resonance the inner eccentricity e\ initially 
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grows linearly at approximately the same rate as it would without the general 
relativistic perturbation. However, the period of the eccentricity oscillation can 
be many times the period of the classical eccentricity oscillations. Although this 
allows the eccentricity to grow to a larger value, the timescale for this growth 
can then be longer than the expected lifetime of the triple in the cluster core. 
For example, with ni2 = 0.006 M we find that the inner binary reaches an 
eccentricity of 0.025 after about 1.5 times the expected lifetime of the triple in 
the core of M4. Thus, even if the system is near resonance, it must probably still 
be residing somewhat outside the core for the secular eccentricity perturbation 
to have enough time to grow to the currently observed value. 

4. Summary 

Our theoretical analysis of the latest timing data for PSR B1620— 26 clearly 
confirms the triple nature of the system. Indeed, the values of all five measured 
pulse frequency derivatives are consistent with our basic interpretation of a bi- 
nary pulsar perturbed by the gravitational influence of a more distant object 
on a bound Keplerian orbit. The results of our Monte-Carlo simulations based 
on the four well-measured frequency derivatives and preliminary measurements 
of short-term orbital perturbation effects in the triple are consistent with the 
complete solution obtained when we include the preliminary measurement of 
the fifth frequency derivative. This complete solution corresponds to a second 
companion of mass rri2 sin i2~7x 10~ 3 M Q in an orbit of eccentricity e2 ~ 0.45 
and semimajor axis 02 ~ 60 AU (orbital period P2 ~ 300 yr). Although the 
present formal la error on is large, we do not expect this solution to change 
significantly as more timing data become available. 

It is possible that the dynamical interaction that formed the triple also 
perturbed the eccentricity of the binary pulsar to the anomalously large value of 
0.025 observed today. However, we have shown that, through a subtle interaction 
between the general relativistic corrections to the binary pulsar's orbit and the 
Newtonian gravitational perturbation of the planet, this eccentricity could also 
have been induced by long-term secular perturbations in the triple after its 
formation. The interaction arises from the near equality between the general 
relativistic precession period of the inner orbit and the period of the Newtonian 
secular perturbations for a low-inclination system. It allows the eccentricity 
to slowly build up to the presently observed value, on a timescale that can be 
comparable to the lifetime of the triple in M4. 

All dynamical formation scenarios have to confront the problem that the 
lifetimes of both the current triple and its parent star-planet system are quite 
short, typically ~ 10 7 — 10 8 yr as they approach the cluster core, where the 
interaction is most likely to occur. Therefore, the detection of a planet in orbit 
around the PSR B1620— 26 binary clearly suggests that large numbers of these 
wide star-planet systems must exist in globular clusters, since most of them will 
be destroyed before (or soon after) entering the core, and most planets will not be 
able to survive long in a wide orbit around any millisecond pulsar system (where 
they may become detectable through high-precision pulsar timing). Although 
a star-planet separation ~ 50 AU may seem quite large when compared to 
the orbital radii of all recently detected extrasolar planets (which are all smaller 
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than a few AU; see Marcy & Butler 1998), one must remember that the current 
Doppler searches are most sensitive to planets in short-period orbits, and that 
they could never have detected a low-mass companion with an orbital period 
3> 10 yr. Similarly, the recent HST search for planetary transits in 47 Tuc by 
Gilliland et al. (2000), which did not detect any planet, was only sensitive to 
very short orbital periods <J 10 d. In addition, it is of course possible that the 
parent system may have been a primordial binary star with a low-mass, brown 
dwarf component, rather than a main-sequence star with planets. 
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